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Abstract. We provide a status report on our current understanding of the mass scales for Pop III. 1 
and Pop III. 2 stars. Since the last review [14], substantial progress has been made both numerically 
and analytically on the late stages of protostellar cloud core collapse, protostar formation and 
accretion, and stellar evolution taking into account cloud core properties and radiative feedback 
effects. Based on this, there are growing indications that primordial stars forming from purely 
cosmological initial conditions (Pop III. 1) were substantially more massive than stars forming in 
preionized gas (Pop III. 2) where HD cooling is important. Different stellar endpoints are predicted 
for these two types of Pop III stars with different chemical enrichment signatures: the former die 
as pair instability supernovae or intermediate mass black holes, whereas the latter die as iron core- 
collapse supernovae, leaving behind neutron star and stellar black hole remnants. We review recent 
simulations which show evidence for binary fragmentation at high densities, and comment on the 
significance of these results. We then summarize an attempt to directly calculate the Pop III. 1 IMF 
taking into account the latest numerical and analytical models. We conclude with suggestions for the 
kind of simulations needed next to continue improving our understanding of Pop III star formation, 
which is a necessary input to understanding high redshift galaxy formation. 
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INTRODUCTION 

The standard model of Pop III star formation holds that they are born massive (~ 
1OOM ) and form as isolated, single stars in the centers of dark matter minihalos of 
mass ~ 1O 6 M at redshifts z=30-15 ([14] and references therein.) The mass scale for 
Pop III stars is fundamentally linked to the cooling properties of primordial gas that viri- 
alizes to T ~ 1000K in such halos. Trace amounts of H2 form via the H channel which 
cools the densest gas to ~ 200K — a temperature set by the excitation temperature of the 
lowest ro vibrational transition. At densities of ~ 10 4 cm~ 3 the H2 level populations of 
assume Boltzmann statistics and cooling becomes independent of density. The Bonner- 
Ebert mass for gas of this density and temperature is about 500M Q . Once H2 cooling has 
accumulated this much mass at these conditions in the central part of the halo, gravita- 
tional collapse ensues. However, due to the low fractional abundance of H2 molecules 
(~ 10~ 4 ), collapse is a quasistatic cooling flow until densities reach about 10 8 cm -3 . At 
this point, gas becomes fully molecular via the 3 -body formation process and the col- 
lapse becomes dynamic. What happens next is fundamental to determining the Pop III 
mass scale and IMF. Does the gas fragment into smaller objects or accrete onto a single 
massive protostar? What are their characteristic masses and mass distributions? 

In this contribution I will review work carried out since my 2008 review [14] that 
examines the high density evolution of the primordial protostellar cloud core, the for- 



mation and growth of the protostar, and stellar evolution through to its endpoint. Recent 
numerical simulations suggest that massive single Pop III stars and also binary stars 
are plausible outcomes. Scenarios are also discussed which, depending on additional 
physics assumptions, may produce more than two fragments in a bound system which 
opens the possibility of dynamic ejection of a low mass primordial star that never ac- 
cretes more mass. The role of radiative feedback on setting the final stellar mass is briefly 
reviewed, but see also Tan et al. (these proceedings). 

FORMATION OF GRAVITATIONALLY COLLAPSING CORES 

We adhere to the nomenclature proposed by McKee at First Stars III [20]. Pop III. 1 stars 
are those that form from pristine initial conditions set by cosmology and primordial gas 
physics. Pop III. 2 stars form from zero metallicity gas that has been pre-processed by 
radiative or mechanical feedback effects, but not by chemical enrichment. Examples of 
the latter would be stars forming in relic HII regions [17], in the presence of a strongly 
photodissociating UV backround [10, 28, 19], or in shocked primordial gas shells [11]. 
Technically, stars forming from gas of such low metallicity that their formation dynamics 
are indistinguishable from Pop III are also included in this class. 



Pop III.1 

A growing list of numerical simulations starting from cosmological initial confirms 
the basic picture outlined in the Introduction. As reviewed in [14], both AMR and SPH 
simulations give results that agree well with one another for the same cosmology, input 
physics, and numerical resolution up to central densities of 10 12 cm~ 3 where they can be 
compared. Evolution to higher central densities is still being explored and is discussed 
below. 

The formation timescale and dynamics for cores forming in the presence of a pho- 
todissociating UVB is not significantly different from the zero UVB case for Jtw < 
10~ 23 erg/cm 2 /s/Hz/ster = 10 _2 /2i [19], and therefore operationally we can include 
them in Pop DLL The evolution of the Lyman- Werner (LW) background due to Pop 
III. 1 sources predict that this level is reached at z ~ 30 [29, 27], and thus we can take 
this as the nominal transition from Pop III. 1 to Pop III. 2. Note that the transition redshift 
is sensitive to Os and the mass of the first stars (through their luminosities). Yoshida 
et al. [29] and Wise & Abel [27] both assumed WMAP1 parameters (<j 8 = 0.9), which 
would overestimate the transition redshift based on the WMAP7 value of Cg = 0.82. It 
is also important to note that proximity effects (i.e. the suppression of H2 formation and 
cooling in clustered minihalos) are ignored in these estimates. The highly clustered na- 
ture of high redshift structure formation potentially makes this an important effect which 
would increase the transition redshift. A lower limit to the transition redshift comes from 
the buildup of the LW background due to normal stars (Pop II) forming in high redshift 
dwarf galaxies. Ahn et al. [3] simulated this process in a 50 Mpc box by considering 
radiation escaping from galaxies with dynamical mass > 1O 8 M (atomic line coolers). 
This simulation explicitly ignores the contribution of minihalos to the LW background. 



They found that the volume averaged /21 reaches 10~ 2 at z ~ 20, which we can take as 
an approximate lower limit for Pop III. 1 formation given the definition above. 

In the absence of a LW background, O'Shea and Norman [18] explored the effect 
of formation redshift on the properties of the collapsing cloud core by simulating 12 
cosmological realizations using the AMR code Enzo. They found that halos virializing 
at higher redshifts (z > 30) systematically produce lower mass cores than z ~ 20 halos 
because the higher virial temperature for a given halo mass translates into more H2 
production and hence more cooling. More cooling leads to lower temperatures and 
accretion rates in the collapsing core, suggesting that Pop III. 1 stars forming at high 
redshift may be less massive than those forming later. Tan et al. (these proceedings) has 
investigated the effect this might have on final stellar masses by applying the accretion 
shutoff model of McKee & Tan [12] (hereafter MT08). They confirm the existence of 
the effect but find that is it weak and dominated by intrinsic scatter. Implications of this 
for the Pop III. 1 IMF are discussed below. 



Pop III.2 

Primordial stars forming from initial conditions modified by radiative and/or mechan- 
ical feedback are called Pop III. 2 stars. Two cases are particularly important: 1) Pop III 
stars forming in the presence of a strong LW background /21 > 10~ 2 (z < 30) from gas 
that has not been preionized or otherwise disturbed, and 2) Pop III stars forming in pre- 
ionized gas that has subsequently recombined and formed H2 and HD molecules, such 
as in a relic HII region. 

In the first case, we imagine a minihalo of characteristic mass 1O 6 M forming signif- 
icantly later and far away from than those forming from the rare (> 3 a) peaks which 
formed Pop III. 1 stars. Yoshida et al. [29] showed that H2 production in minihalos is 
suppressed for /21 > 10 2 . Using high resolution AMR simulations O'Shea & Nor- 
man [19] (hereafter ON08) showed that collapsing core formation is not eliminated for 
J21 > 10~ 2 , but only substantially delayed. They showed that even for values as large as 
721 = 1 core collapse eventually happens because although the abundance of H2 is very 
low at these background levels, the cooling rate per molecule is a strongly increasing 
function of temperature. At some point the halo becomes massive enough through ac- 
cretion and mergers that its virial temperature becomes high enough so that the higher 
cooling rate per particle multiplied by a long evolutionary time (~ 200 Myr) removes 
enough thermal energy to trigger gravitational collapse. Such "Lyman- Werner delayed 
halos" accumulate more mass in the core and of higher temperature (see OSShea & 
Whalen, these proceedings) than Pop III. 1 star-forming halos. Once collapse begins in a 
delayed halo, the accretion rates are higher as well, suggesting that either the Pop III. 2 
star will be more massive than Pop III. 1 stars if only a single star is formed, or that 
several Pop III.2 stars will form. Calculations of this sort have not been carried to a high 
enough density to know which of these two possibilities is more likely. 

The second case — primordial stars forming in recombining gas — is expected to be 
common, and possibly even the dominant mode of Pop III star formation. The reasons 
for this are several. First, due to the clustering in CDM models, the second halo to form 



a Pop III star is very likely to be near the first one to do so in a given volume of the 
universe. This is borne out by simulations [17, 13, 31]. Second, Pop III stellar lifetimes 
are only about 1% of the Hubble time during the Pop III epoch, and recombination times 
are also short compared to the Hubble time for z > 10, meaning HII regions come and 
go quickly in a given region of space, even as the ionized volume fraction slowly grows. 
Third, Pop III star formation rate density, and the ionized volume fraction secularly 
increases toward lower redshifts, meaning that the likelihood a Pop III star forms in this 
way also increases toward lower redshifts. 

Primoridal stars forming in relic HII regions form from initial conditions that are 
thermally, chemically, and structurally altered relative to the Pop III. 1 case. The most 
careful simulations of this scenario were carried out by Mesinger et al. [13] and Yoshida 
et al. [31]. They found that the high electron abundance in the recombining gas catalyzed 
rapid formation of H2 molecules which roughly compensated the Jeans smoothing of 
the gas during the ionized phase. Ignoring the effects of HD formation and cooling, they 
found that core properties and accretion rates were very similar to the Pop III. 1 case. 
However, when HD cooling was included, they found lower core temperatures (100K 
vs. 200K) which lowered the core accretion rates substantially. Simple accretion time 
arguments suggested typical final stellar masses of 40 — 6OA/ , rather than 100 — 3OOM 
characterizing Pop III. 1 stars. 

FROM COLLAPSING CORES TO PROTOSTARS 

Additional physical processes must be included to evolve the collapsing cloud core to 
stellar density. These include: 1) heat of formation of Hi and change of the adiabatic 
index above n ~ 10 8 cm~ 3 ; 2) radiative transfer of the H2 cooling lines above n ~ 
10 12 cm~ 3 ; 3) radiative cooling due to collisionally induced emission (CIE) above n ~ 
10 14 cm~ 3 including continuum opacity effects; and 4) cooling due to H% dissociation 
above n ~ 10 15 cm~ 3 . Omukai & Nishi (1998)[16] were the first to include all of these 
processes in ID spherically symmetric collapses. A decade later, Yoshida and colleagues 
achieved it in 3D [32] starting from proper cosmological initiation conditions-truly 
a milestone in the subject. Radiative transfer effects were modeled using an escape 
probability formalism calibrated against full radiative transfer calculations in ID. 

Fig. 1 shows structures on four different lengthscales from the final snapshot of 
their simulation which spans a range of scales of 13 orders of magnitude! The newly 
elucidated structures in this simulation are the centrifugally supported, Toomre unstable 
disk in Fig. lc, and the hydrostatic central protostar in Fig. Id. The latter is formed with 
a mass of 1O~ 2 M and is accreting at a rate between 0.01 and 0.1 M Q /yr, implying 
that it will grow to a mass of tens of solar masses in just 10 4 yr. The simulation, which 
resolves the Jeans mass everywhere, shows no evidence for fragmentation. 

Time-explicit numerical simulations such as this cannot follow the accretion phase 
because of the very small timestep imposed by the Courant condition. However, the ac- 
cretion rate history onto the central protostar can be estimated from a snapshot of the 
mass density and velocity distributions in the envelope, and then fed into a ID protostel- 
lar evolution code. Ignoring late stage fragmentation, Yoshida et al. [31] estimate a Pop 
111.1(111.2) enters the ZAMS with a mass of ~ 1OO(4O)M , respectively. 
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FIGURE 1. A zoom-in showing structures on diverse length scales from the first 3D hydrodynamic 
cosmological simulation to form a hydrostatic, accreting primordial protostar. From [32]. 



It is important to note that Pop III stars continue to accrete after the onset of hydrogen 
burning. Their final mass will be set either by the mass available to accrete, or by 
radiative feedback effects [12]. Interestingly, in the case of a Pop III. 2 star forming in a 
relic HII region, [31] found that the cloud core was only about 4OM , suggesting that 
this number may be close to the final stellar mass. Pop III. 1 stars, and Pop III. 2 stars 
forming in Lyman-Werner delayed halos are surrounded by much more massive cloud 
cores and may accrete until radiation terminates the accretion. This is discussed further 
below. 

FRAGMENTATION 

Core fragmentation due to turbulence 

The collapsing cores of primordial minihalos are characterized by a state of subsonic 
turbulence, and therefore there is a ready source of perturbations to trigger fragmenta- 
tion if the gas is thermally or chemo-thermally unstable [22]. Fragmentation was not 
observed in the cosmological simulations of Abel et al. [1, 2] and Yoshida et al. [30, 32], 
although it was observed in the idealized collapse simulations of Bromm et al. [5] and 
Clark et al. [7]. Clearly, initial conditions matter. We still have rather small samples 
of self-consistent, high dynamic range simulations of Pop III star formation which in- 
form our standard model. O'Shea & Norman [18] simulated four cosmic realizations in 
each of three different box sizes and found no fragmentation up to central densities of 
10 12 cm~ 3 . However, Turk et al. [24] found one case out of 5 higher density evolutions 



in which the core fragments into two subcores, each of which were in a state of freefall 
collapse. With a separation of nearly 1000 AU, Turk et al. suggest the final outcome is a 
long period Pop III binary. Turk (these proceedings) reports a 20% binary incidence rate 
in a larger sample of 25 simulations. 



Late fragmentation due to disk instabilities 

The formation of a rotationally supported, fully molecular disk as in Fig. lc creates an 
environment for fragmentation. The disk will fragment if it becomes Toomre unstable. 
Toomre instability depends on the gas surface density k through the Q parameter: 
Q = c s k/ kGL where K is the epicyclic frequency. If material is added to the disk faster 
than it can be deposited onto the central object, the surface density will grow to the point 
where the disk becomes gravitationally unstable. The strong m = 2 mode seen in Fig. 
lc is the first manifestation of global Toomre instability. If Q drops substantially below 
unity, the disk becomes locally gravitationally unstable and can fragment into multiple 
objects. 

A technical impediment to exploring possible disk fragmentation in a numerical 
simulation like that shown in Fig. 1 is the telescoping down of the timestep as the central 
object becomes denser and hotter. One approach to avoiding the "Courant catastrophe" 
is to replace the central object with a sink particle, thus regularizing the simulation. Sink 
particles were introduced by Bate et al. [4] in the context of Galactic star formation 
simulations, and first employed for primordial star formation by [6]. At an operational 
level, all gas above a given density is converted into a sink particle thus placing a floor 
on the timestep. The sink gains mass from its environment as gas particles cross its 
gravitational radius. 

Stacy et al. [23] employed sink particles to study the evolution of the accretion disk 
surrounding a Pop III. 1 protostar in a standard cosmological setup. They chose as their 
threshold density for creating a sink 10 12 cm~ 3 by which time the gas is fully molecular 
but before the additional physical processes mentioned above become important. They 
integrated for 5000 yr after sink creation, during which time the sink grew in mass from 
0.7 to 43 Mq. They terminated the calculation before radiative feedback effects become 
important. They found that a centrifugally supported disk forms around the protostar 
which grows to a size of ~ 2000 AU and mass ~ 4OM in 5000 yr. The disk fragments 
into 5 objects, two of which are massive and long-lived (the central star and one at 700 
AU). 

One needs to be careful accepting this result at face value. The first concern is the 
sensitivity of the results to the threshold density for making a sink. A higher/lower 
density will likely result in a smaller/larger disk (this effect has been seen in AMR 
simulations, where an artificial pressure is used to regularize the simulation.) Assuming 
the disks fragment in a similar way, the separation of the multiple protostars would be 
sensitive to the threshold density. This is borne out by the higher resolution simulations 
presented by Clark (these proceedings) which used a threshold density of 10 15 cm~ 3 . The 
second concern is the correctness of the dynamical evolution of the sink particles once 
formed. While a sink particle treatment for the central protostar seems "safe" because 
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FIGURE 2. Left: Schematic of the evolutionary phases of an accreting Pop III protostar. Right: Proto- 
stellar accretion rate versus accreted mass including radiative feedback. When the disk photoevaporation 
rate equals the accretion rate, accretion stops. From [12]. 

it sits in the center of the potential well and accretes rapidly, condensations forming 
in the disk are subject to very different stresses (gravitational torques, hydrodynamic 
shear). Converting extended condensations into sinks gives them an indestructability 
that may not be realistic. Moreover, whether the condensations migrate inward and 
accrete onto the central star, or migrate outward to become a stable binary/multiple 
system depends on modeling the torques and hydrodynamic drag accurately. Despite 
these caveats, however, these simulations do lend strong support to the notion that at 
least some and perhaps all Pop III stars form as binaries. 

FROM PROTOSTARS TO STARS TO STELLAR ENDPOINTS 
Halting accretion: radiative feedback 

Because Pop III stars are predicted to be ultraluminous [21], it is important to consider 
the effects of radiation on the accreting envelope. McKee & Tan [12] studied this within 
the context of analytic models of Pop III star formation. The model assumes that the Pop 
III protostar accretes from a disk whose properties are estimated from the conditions 
found in the collapsing cloud cores from the cosmological simulations of Abel et al. 
[2] (hereafter ABN02). They considered various forms of radiative feedback, including 
H2 photodissociation, Lyman a and Lyman continuum pressure, HII region breakout, 
and accretion disk photoevaporation. They found that the first three of these effects are 
not significant at any time, and the last two only become important at later times as the 
protostar becomes more massive and luminous. 

Fig. 2 illustrates the basic mechanism that shuts off disk accretion. As the accretion 
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FIGURE 3. Evolution of Pop III protostars and stars to their endpoints. Left: without radiative feedback; 
Right: with radiative feedback. Cases Y-l to Y-5 are Pop III. 1 protostars with accretion rates reduced by 
varying amounts due to rotation. Case YII is a Pop III. 2 protostar. Cases M-l to M-3 are Pop III. 1 stars 
varying the ionized gas temperature and accretion rate. From [15]. 



rate from the envelope diminishes, and the star's UV luminosity rises, an HII region 
breaks out above and below the disk. The expanding HII region blocks accretion above 
and below the disk, reducing the accretion rate below the no feedback value, and bathes 
the disk in ionizing radiation from recombinations within the nebula. The disk begins to 
photoevaporate. When the disk photoevaporation rate equals the accretion rate, inflow is 
halted and the star attains its final mass. The final mass depends on the entropy of the gas 
accreting from large radius, as well as its specific angular momentum. Fig. 2b shows the 
fiducial case from MT08 taking values from the Pop III. 1 simulation of ABN02. They 
predict a final mass of ~ 14OA/ . 



Pop III stellar evolution 

Ohkuba et al. [15] used accretion histories from the cosmological simulations of 
[31, 32] to carry out Pop III. 1 and Pop III.2 stellar evolution calculations through to their 
end points. To consider some variation in the parent cloud, they parameterize angular 
momentum and the radiative feedback effects of MT08. Ignoring radiative feedback, 
they find Pop III. 1 stars die as very massive stars (VMS) in the range 300 — 1000M Q , 
depending on the angular momentum of the cloud. These stars produce intermediate 
mass black hole (IMBH) remnants and no chemical enrichment. Note this mass range 
does not include the range 140 — 260M Q predicted to produce pair instability supernovae 
(PISN)[8]. However, when radiative feedback is included, they find that Pop III. 1 stars 
die with masses in the range 60 — 320M Q , depending on angular momentum of the 
cloud. This mass range includes Type II supernovae, PISN, and IMBHs as possible 
outcomes. Clearly, radiative feedback effects are fundamental in determining the Pop III 
IMF and chemical signatures. Finally, they find that Pop III.2 stars die as core-collapse 
supernovae in the mass range 40 — 6OM , consistent with earlier estimates [13, 31]. 
These stars produce stellar black holes and some chemical enrichment. 



Given the findings described above, what can we say about the Pop III IMF? Ignoring 
binarity for the moment, it is now quite clear that Pop III. 1 stars are considerably more 
massive than Pop III. 2 stars forming in relic HII regions. Table 1 in MT08 quote values 
of the photoevaporation limited mass in the range of ~ 60 — 32OM varying parameters 
over a reasonable range, with a fiducial mass of 14OA/ . Ohkubo et al. [15] find similar 
values for the stellar endpoint masses, although this is not surprising since they used 
fits to MT08 rates as input to their calculations. Tan et al. (these proceedings) used the 
ensemble of 12 simulations of [18] to estimate the Pop III. 1 IMF. They confirmed that 
stars forming at lower redshift are more massive, but caution that the trend is weak and 
there is a large scatter. The key factor setting the final mass is the entropy parameter 
for the gas in the Jeans unstable core. High masses are obtained in higher temperature 
cores. They predict final masses in the range 80 — 6OOM , with a mean(median) value 
of 25O(215)M , respectively. Based on this small sample, they compute a skewed IMF 
that has a peak at ~ M , falling off very rapidly to lower masses, and more gradually to 
higher masses. 

Although the sample of cosmological simulations of Pop III.2 stars forming in relic 
HII regions is even smaller [17, 13, 31], there is a consensus that a lower value of 
40 — 6OM is more appropriate. This mass reflects the smaller mass of the cloud core 
to accrete from, and the lower accretion rates that results from HD cooling the gas to 
< 100K. At these mass scales, radiation feedback is less important in setting the mass 
scale. 

Pop III. 2 stars forming in Lyman- Werner delayed halos should fall in the mass range 
of the Pop III. 1 stars mentioned above, and possibly extend it to higher masses [19]. The 
reason is that delay translates into higher virial masses and temperatures for the minihalo 
before gravitational instability sets in. Once it sets in, accretion rates are higher. ON08 
observe a discontinuous jump in the virial mass/temperature of the minihalo at the onset 
of gravitational instability between J LW = 10~ L5 /2i and 10 V21 (see Fig. 1 of O'Shea, 
these proceedings, or Fig. 3 of ON08.) These values are (2.5 x 1O 6 M ,4OOOK) and 
(1O 7 M , 9000K), respectively. The mass of the collapsing cores when the central density 
reaches 10 10 cm~ 3 are 3 x 1O 3 M and 6 x 1O 3 M , respectively. Provided such objects can 
form in the universe, we might expect a discontinuous change in the characteristic stellar 
masses that are formed at these values of Jlw- Once collapse begins, the H2 fraction at 
the onset of collapse becomes unimportant, only the mass and entropy of the collapsing 
core. Although no fragmentation was seen in the simulations to a central density of 
10 10 cm~ 3 , we have seen that fragmentation can occur at higher densities. With the larger 
supply of infalling material, it is possible that the outcome is a small cluster of ~ 1OOM 
stars. 

The possibility of binary formation does not change the basic picture outlined above, 
but merely divides the mass scales by two. If final masses are determined by the amount 
of gas available to accrete onto the binary protostars, as in relic HII regions, then masses 
of 20 — 3OM would be predicted. However, Pop III. 1 stars and Pop III.2 stars forming 
in Lyman- Werner delayed halos have a larger reservoir of gas to accrete from. Final 
masses may be set by a complex interplay between radiative feedback and high density 
fragmentation effects. More gas may simply mean more stars of the characteristic mass 
~ 1OOM set by photoevaporation limited disk accretion. 



WHAT SIMULATIONS ARE NEEDED NOW? 



We need larger samples of cores from cosmological initial conditions reaching central 
densities of at least 10 12 cm~ 3 in various Lyman- Werner fluxes to explore the turbulent 
state of the gas and core multiplicities, as well as to quantify the entropy and rotational 
parameters entering in the MT08 model. From this sample an improved estimate of the 
Pop III. 1 IMF can be obtained using the procedures of Tan et al. (these proceedings). 

We need a statistical sample of cores from cosmological initial conditions forming in 
relic HII regions, taking the detailed minihalo photoevaporation dynamics and Lyman- 
Werner background into account. Idealized studies of minihalo photoevaporation in- 
cluding multifrequency radiative transfer and nonequilibrium primordial gas chemistry 
by Whalen et al. [25, 26] show a complicated set of outcomes depending upon the Pop 
III. 1 star's luminosity and spectrum, and the Pop III.2 halo's distance and evolutionary 
state (O'Shea & Whalen, these proceedings). Competing effects (e.g., H2 photodissocia- 
tion, X-ray preionization and H2 formation, and hydrodynamic cloud crushing) roughly 
cancel out in many cases, leaving neighboring minihalos is a state of collapse after the 
I-front has swept over it. Given the expectation that most of Pop III stars from in the 
vicinity of other Pop III stars and other ionizing sources, this is a rich class of models 
that need further investigation. 

Sink particles are clearly a valuable tool in investigating the late time, high density 
evolution of the central region of the cloud. A careful study of disk formation and 
fragmentation varying resolution and sink recipes would be invaluable to improving 
our confidence in the results, and elucidating the question of stellar multiplicity. The 
simulations by Clark et al. (these proceedings) are a step in the right direction. 

3D radiation hydrodynamical simulations of MT08 model coupled to a radiating 
protostar evolution model modeled as a sink particle would be an important next step 
in understanding the details of accretion shutoff. The simulations would be analogous to 
the massive star formation simulations of Krumholz et al. [9], but include the transport 
of ionizing radiation. These simulations could profitably be done in idealized setups or 
in fully cosmological settings. 

Finally, coupled models of the evolution of the Lyman-Werner background and hy- 
drogen reionization in cosmological volumes with adequate resolution to resolve the 
minihalo population would be numerically challenging but extremely valuable to assess 
the duration of the Pop III. 1 and Pop III.2 epoch, and the impact of primordial stars of 
all types on the initial conditions for galaxy formation. 
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